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Compress Images 


Mathematics is providing a novel technique for achieving 
compression ratios of 10,000 to 1—and higher 
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THE NATURAL WORLD is filled with 
intricate detail. Consider the geometry 
on the back of your hand: the pores, the 
fine lines, and the color variations. A 
camera can capture that detail and, at 
your leisure, you can study the photo to 
see things you never noticed before. Can 
personal computers be made to carry out 
similar functions of image storage and 
analysis? If so, then image compression 
will certainly play a central role. 

The reason is that digitized images— 
images converted into bits for processing 
by a computer—demand large amounts of 
computer memory. For example, a high- 
detail gray-scale aerial photograph might 
be blown up to a 3%2-foot square and then 
resolved to 300 by 300 pixels per square 
inch with 8 significant bits per pixel. 
Digitization at this level requires 130 
megabytes of computer memory—too 
much for personal computers to handle. 

For real-world images such as the aeri- 
al photo, current compression techniques 
can achieve ratios of between 2 to 1 and 
10 to 1. By these methods, our photo 
would still require between 65 and 13 
megabytes. 

In this article, we describe some of the 
main ideas behind a new method for 
image compression using fractals. The 
method has yielded compression ratios in 
excess of 10,000 to 1 (bringing our aerial 
photo down to a manageable 13,000 
bytes): The color pictures in figures 1 
through 5 were encoded using the new 
technique; actual storage requirements 
for these images range from 100 to 2000 
bytes. 


A mathematics research team at the 


Georgia Institute of Technology is devel- 
oping the system, with funding provided 
by the Defense Advanced Research Proj- 
ects Agency (DARPA) and the Georgia 
Tech Research Corporation (GTRC). Our 
description is necessarily simplified, but 
it will show you how a fractal image- 
compression scheme operates and how to 
use it to create exciting images. 


Describing Natural Objects 
Traditional computer graphics encodes 
images in terms of simple geometrical 
shapes: points, line segments, boxes, cir- 
cles, and so on. More advanced systems 
use three-dimensional elements, such as 
spheres and cubes, and add color and 
shading to the description. 

Graphics systems founded on tradi- 
tional geometry are great for creating pic- 
tures of man-made objects, such as 
bricks, wheels, roads, buildings, and 
cogs. However, they don’t work well at 
all when the problem is to encode a sun- 
set, a tree, a lump of mud, or the intricate 
structure of a black spleenwort fern. 
Think about using a standard graphics 
system to encode a digitized picture of a 
cloud: You'd have to tell the computer the 
address and color attribute of each point 
in the cloud. But that’s exactly what an 
uncompressed digitized image is—a long 
list of addresses and attributes. 

To escape this difficulty, we need a 
richer library of geometrical shapes. 
These shapes need to be flexible and con- 
trollable so that they can be made to con- 
form to clouds, mosses, feathers, leaves, 
and faces, not to mention waving sun- 
flowers and glaring arctic wolves. Fractal 


geometry provides just such a collection 
of shapes. For a hint of this, glance at the 
pictures in The Fractal Geometry of Na- 
ture by Benoit Mandelbrot, who coined 
the term fractal to describe objects that 
are very “fractured” (see references for ~ 
additional books and articles). Some ele- 
mentary fractal images accompany this 
article. 

Using fractals to simulate landscapes 
and other natural effects is not new; it has 
been a primary practical application. For 
instance, through experimentation, you 
find that a certain fractal generates a pat- 
tern similar to tree bark. Later, when you 
want to render a tree, you put the tree- 
bark fractal to work. : 

What is new is the ability to start with 
an actual image and find the fractals that 
will imitate it to any desired degree of ac- 
curacy. Since our method includes a com- 
pact way of representing these fractals, 
we end up with a highly compressed data 
set for reconstructing the original image. 


Overview of Fractal Compression 
We start with a digitized image. Using 
image-processing techniques such as 
color separation, edge detection, spec- 
trum analysis, and texture-variation 
analysis, we break up the image into seg- 
ments. (Some of the same techniques 
continued 
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Figure 1: /FS-encoded color image of 
three-dimensional ferns (4 
transformations, 100 bytes). 





Figure 2; [FS-encoded color photo of 
Black Forest, color set adjusted to give 
winter tones (120 transformations, 
2000 bytes). 





T 


Figure 3: [FS-encoded color photo of a 
Bolivian girl (120 transformations, 2000 
bytes). 
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form the basis for the automatic coloring 
of black-and-white motion pictures.) A 
segment might be a fern, a leaf, a cloud, 
ora fence post. A segment can also be a 
more complex collection of pixels: A sea- 
scape, for example, may include spray, 
rock, and mist. 

We then look up these segments in a li- 
brary of fractals. The library doesn’t 
contain literal fractals; that would require 
astronomical amounts of storage. In- 
stead, our library contains relatively 
compact sets of numbers, called iterated 
function system (IFS) codes, that will re- 
produce the corresponding fractals. Fur- 
thermore, the library’s cataloging system 
is such that images that look alike are 
close together: Nearby codes correspond 
to nearby fractals. This makes it feasible 
to set up automated procedures for 
searching the library to find fractals that 
approximate a given target image. A 
mathematical result known as the Collage 


~ Theorem (more on that later) guarantees 


that we can always find a suitable IFS 
code—and gives a method for doing so. 

Once we have looked up all the seg- 
ments in our library and found their IFS 
codes, we can throw away the original 
digitized image and keep the codes, 
achieving our compression ratio of 
10,000 to 1—or even higher. 


Iterated Function Systems 
We start by explaining how a set of IFS 
codes can approximate a natural image. 
IFS theory is an extension of classical 
geometry. It uses affine transformations, 
explained below, to express relations be- 
tween parts of an image. Using only these 
relations, it defines and conveys intricate 
pictures. With IFS theory, we can de- 
scribe a cloud as clearly as an architect 
can describe a house. 
By studying the following sections, 





= 


Figure 4: [FS-encoded color photo of 
the Monterey coast (60 transformations, 
100 bytes). 





you should be able to encode and decode 
fascinating black-and-white image seg- 
ments, such as leaf skeletons, tree 
shadows, spirals, and thunderheads. You 
should also obtain an overview of how a 
fully automated fractal compression sys- 
tem operates. 

Affine transformations can be de- 
scribed as combinations of rotations, 
scalings, and translations of the coordi- 
nate axes in n-dimensional space. An ex- 
ample in two dimensions is 


Wx,y) = (Axt My+1, 4x+"%y+2), 


which can also be written in matrix form 
as 


"Ul [3s 2] B] Led. 


This transformation moves the point (0,0) 
to (1,2) and moves (—1,0.5) to (0.625, 
2). To confirm your understanding of the 
idea, you should work out where it moves 
the point (1,1). We denote this transfor- 
mation by W; the notation W(S) denotes 
the subimage of W on a set of points S. 

Now let’s see what W does to a picture 
of a smiling face, F, lying on the x,y plane 
(see figure 6). The result is a new, 
squeezed face W(F). The affine transfor- 
mation has deformed and moved the face. 
Notice that the eyes in the transformed 
face W(F) are closer together than they 
are in F. We say that the transformation W 
is contractive: It always moves points 
closer together. 

Another example of a contractive af- 
fine transformation is shown in figure 7. 
This time it acts on a leaf to produce a 
new, smaller leaf. 

The general form for an affine trans- 
formation is 





Figure 5: /FS-encoded color image 
from A Cloud Study (30 
transformations, 500 bytes). 


Photographs © 1987 Georgia Tech Research Corp. Figure 1: Barnsley; Figures 2 and 3: Barnsley, 
Jacquin, Malassenet, Reuter, Sloan; Figures 4 and 5: Barnsley, Jacquin, Reuter, Sloan. 
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“5]-[22] 5)+[F] 
co 


where the coefficients a, b, c, d, e, and f 
are real numbers. 

If we know in advance the translations, 
rotations, and scalings that combine to 
produce W, we can generate coefficient 
values as follows: 


a=rcos0,b = —ssing, 
c=rsin6,d=scos¢, 


where r is the scaling factor on x, s is the 
scaling factor on y, 6 is the angle of rota- 
tion on x, ¢ is the angle of rotation on y, e 
is the translation on x, and f is the transla- 
tion on y. 

How can you find an affine transfor- 
mation that produces a desired effect? 
Let’s show how to find the affine trans- 
formation that takes the big leaf to the lit- 
tle leaf in figure 7. We wish to find the 
numbers a, b, c, d, e, and f for which the 
transformation W has the property 


W (big leaf) = little leaf. 


Begin by introducing x and y coordinate 
axes, as already shown in the figure. 
Mark three points on the big leaf (we’ve 
chosen the leaf tip, a side spike, and the 
point where the stem joins the leaf) and 
determine their coordinates (a;,a2), 
(81,82), and (1,72). Mark the corre- 
sponding points on the little leaf and 
determine their coordinates (@,,a2), 
(81,82), and (1 »Y2) respectively. 

Determine values for the coefficients 
a, b, and e by solving the three linear 
equations 


Table 1: [FS codes for a Sierpinski triangle. 


b 


b 





M1a+ambt+e=a, (1) 
Biat+Prb+e=8,, (2) 
matybt+e=jn, (3) 


and find c, d, and f in similar fashion 
from these equations: 


ac+tad+f=a, (4) 
Bict+Bd+f=B,, (5) 
metydt+f=. (6) 


We recommend the use of an equation 
solver such as TK Solver Plus (Universal 
Technical Systems, Rockford, Illinois) or 
Eureka (Borland International, Scotts 
Valley, California) for finding the coeffi- 
cient values. Doing it manually can be 
tedious. 

Now that we know what a contractive 





Figure 6: An affine transformation W 
moves the smiling face F to a new face 
W(F). The transformation is called 
contractive because it moves points 





closer together. 
continued 


ed: (a; 2) 


| 
(¥1»¥2) 


Figure 7: Two ivy leaves fix an affine transformation W. 


Table 3: IFS codes for a fern. 
a b c d 


0 0 0 0.16 
0.2 -0.26 0.23 0.22 
-0.15 0.28 0.26 0.24 
0.85 0.04 -0.04 0.85 





Table 4: FS codes for fractal tree. 


WwW a b c d 


0 0 0 0.5 
0.1 0 0 0.1 
0.42 -042 0.42 0.42 
0.42 0.42 -0.42 0.42 
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affine transformation is and how to find 
one that maps a source image onto a de- 
sired target image, we can describe an 
iterated function system. An IFS is a col- 
lection of contractive affine transforma- 
tions. Here’s an example of an IFS of 
three transformations: 


w{5]=[88 93] [5] +[8] 
w.{3]=[38 93] [2] +[8) 
mL T-[85 08] B]+L? 


Each transformation must also have an 
associated probability, p,, determining its 
“importance” relative to the other trans- 
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formations. In the present case we might 
have p,, p2, and p3. Notice that the proba- 
bilities must add up to 1. That is, p; + 
Pit ps=1. 

Of course, the above notation for an 
IFS is cumbersome. Table 1 expresses 
the same information in tabular form. 
Other examples of IFS codes are given 
in tables 2 through 4. Notice that an IFS 
can contain any number of affine 
transformations. 


The Random Iteration Algorithm 
Now let’s see how to decode an arbitrary 
IFS code using the random iteration 
method. Remember that in general an 
IFS can contain any number, say m, of af- 
fine transformations, W,, W2, Ws,..., 
W,,, each with an associated probability. 
The following code summarizes the 
method: 





8: The result of applying the random iteration algorithm to the IFS code in 


table 1. It is called the Sierpinski triangle. 


(i) Initialize: x=0, y=0. 

(ii) For n=1 to 2500, do steps (iii)-(vii). 
(iii) Choose k to be one of the numbers 1, 
2,...,m, with probability p,. 

(iv) Apply the transformation W, to the 
point (x,y) to obtain (x,y). 

(v) Set (x,y) equal to the new point: x=x, 
ye ye 

(vi) Ifn>10, plot (x,y). 

(vii) Loop. 


Applying this procedure to the trans- 
formation in table 1 produces the figure 
shown in figure 8—a fractal known as the 
Sierpinski triangle. Increasing the num- 
ber of iterations n adds points to the 
image. Figure 9 shows the result of the 
random iteration algorithm applied to the 
data in table 3, at several stages during 
the process. By increasing the scale fac- 
tor used in plotting, you can zoom in on 
the image (see figure 10). The text box on 
page 221 contains a BASIC implementa- 
tion of the method with additional com- 
ments on programming. 

You may wonder why the first 10 
points are not plotted (step (vi)). This is 
to give the randomly dancing point time 
to settle down on the image. It is like a 
soccer ball thrown onto a field of expert 
players: Until someone gains control of 
the ball, its motion is unpredictable, or at 
least is independent of the players’ ac- 
tions. But eventually a player gets the 
ball, and its motion then becomes a direct 
result of the skill of the players. The fact 
that our transformation is contractive 
guarantees that the “ball” will eventually 
get to one of the “players,” and that it 
will stay under control after that. 

How do we know that the random iter- 

continued 





Figure 9: A fern appears when the random iteration algorithm is applied to the IFS code in table 3. 
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Circle 244 on Reader Service Card for MS DOS Products. 


(All others: 245) — 


ation algorithm will produce the same 
image over and over again, independent 
of the particular sequence of random 
choices that are made? This remarkable re- 
sult was first suggested by computer-graph- 
ical mathematics experiments and later 
given a rigorous theoretical foundation by 
Georgia Tech mathematician John Elton. 


The Collage Theorem 
Our next goal is to show a systematic 
method for finding the affine transforma- 
tions that will produce an IFS encoding of 
a desired image. This is achieved with the 
help of the Collage Theorem. 

To illustrate the method, we start from 
a picture of a filled-in square S in the x,y 
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plane, with its vertices at (0,0), (1,0), 
(1,1), and (0,1) (see figure 11). The ob- 
jective is to choose a set of contractive af- 
fine transformations, in this case W,, W2, 
W,, W,, so that S is approximated as well 
as possible by the union of the four sub- 
images W,(S)U W,(S)U W,(S) U W,(S). 
Figure 11 shows, on the left, S together 
with four noncovering affine transforma- 
tions of it; on the right, the affine trans- 
formations have been adjusted to make 
the union of the transformed images cover 
up the square. 

To find the coefficients of these trans- 
formations, we use the method described 
earlier in the section on iterated function 
systems, leading to simultaneous equa- 


tions 1 through 3 and 4 through 6. The 
values one finds in the present case are 
given in table 2. When the random itera- 
tion algorithm is applied to this IFS code, 
the square is regenerated. 

The preceding example typifies the 
general situation: You need to find a set 
of affine transformations that shrink dis- 
tances and that cause the target image to 
be approximated by the union of the af- 
fine transformations of the image. The 
Collage Theorem says that the more ac- 
curately the image is described in this 
way, the more accurately the transforma- 
tions provide an IFS encoding of it. 

Figure 12 provides another illustration 
of the Collage Theorem. At the bottom 
left is shown a polygonalized leaf bound- 
ary, together with four affine transforma- 
tions of that boundary. The transformed 
leaves taken together do not form a very 
good approximation of the leaf; in conse- 
quence, the corresponding IFS image 
(bottom right), computed using the ran- 
dom iteration algorithm, does not look 
much like the original leaf image. How- 
ever, as the collage is made more accu- 
rate (upper left), the decoded image 
(upper right) becomes more accurate. 

So, there’s a fundamental stability 
here. You don’t have to get the IFS code 
exactly right in order to capture a good 
likeness of your original image. More- 
over, the IFS code is robust: Small per- 
turbations in the code will not result in 
unacceptable damage to the image. In 
each of the above examples, we have used 
four transformations to encode the image. 
However, any number can be used. 

For example, the spiral image in figure 
13 can be encoded with just two contrac- 
tive affine transformations. See if you can 
find them. Then determine the IFS trans- 
formation coefficients and input them to 
the random iteration algorithm to get the 
spiral back again. 





Figure 11: The collage theorem is used to encode a classical square S. The correct IFS code is obtained when the four affine 
transformations of S cover S, as shown on the right. 
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Assigning Probabilities 

Once you have defined your transforma- 
tions, you need to assign probabilities to 
them. Different choices of probabilities 
do not in general lead to different images, 
but they do affect the rate at which vari- 
ous regions or attributes of the image are 


IFS Decoding 
in BASIC 


t isting A is a BASIC implementation 
of the random iteration algorithm. 
It includes the data for the Sierpinski tri- 
angle, but you can use it to process any 
IFS tables. In particular, you will want 
to try the data in tables 2, 3, and 4. Be 
sure to set the variable m correctly; it 
tells the program how many transforma- 
tions are in the IFS. 

It is also essential that the probabili- 
ties in p() add up to 1. For speed, the 
transformations should be listed in de- 
scending order of probability: the high- 
est probability transformation first, and 
the lowest probability last. 

The program includes variables for 
rescaling and translating the origin to 
accommodate the range of the points be- 
ing plotted to the limits of your screen. If 
the image is too wide, decrease xscale; 
if the points are too close horizontally, 
increase xscale. Adjust yscale simi- 
larly to get a good vertical point spread. 
To move the image, adjust xoffset and 
yoffset. 

You can do these adjustments by trial 
and error: Run the program; interrupt it 
and change the offsets and scale factors; 
and run it again. Or, you can replace the 
plot command pset with a command to 
print the values of x and y and run the 
program to get an exact idea of the range 
of points being plotted, so you can adjust 
the scale and offsets more precisely. 

Another way to arrange the program 
is to have it read all the data—m, a(), 
b{),.,eC),.4(),,e0), £C),.p0), 
xscale, yscale, xoffset, and yoff- 
set—from a disk file specified by the 
user. Instead of reading in the coeffi- 
cients a, b, c, and d, you may want to 
read in angles 6 and ¢ and scale factors r 
and s, and then calculate the 
coefficients. 

The random iteration method is com- 
putation-intensive, so we recommend 
use of a compiler such as Microsoft’s 
QuickBASIC or Borland’s Turbo 
BASIC. If your computer has a floating- 
point coprocessor and your compiler 
supports one, so much the better. 





filled in. Let the affine transformations amount of time that the randomly dancing 
W, corresponding to an image I be point should spend in the subimage VW, is 


approximately equal to 


m()-[eai G)-Ee). Se 


area of I 


where i = 1, 2, 3,..., n. Then the continued 


Listing A: A BASIC program demonstrating the use of the random 
iteration algorithm to reconstruct an IFS-compressed image. 


18 


‘Allow for a maximum of 4 transformations in the IFS 
DIM a(4), b(4), c(4), d(4), e(4), £(4), pl4) 


‘Transformation data, Sierpinski triangle 

‘First comes the number of transformations 

‘then the coefficients a through f and probability pk 
‘The values for pk should be in descending order. 
DATA 3 

DATA .5,8,8, .5,8,8, .34 

DATA .5,8,9,.5,1,8, .33 

DATA. 975010) 6D 2 Decoy 425 


‘Read in the data 

READ m 

pt = 8 ‘Cumulative probability 

FOR j = 170m 
READ a(j), b(j), c(j), d(j). e(5), £(5), pk 
pt = pt + pk 
p{j). = pt 

NEXT j 


‘Set up for Graphics 

SCREEN 3 ‘Select graphics screen 
xscale = 358 ‘Map [8,1] onto [8,358] 
yscale = 325 ‘Map [8,1] onto [8,325] 
xoffset = @ 

yoffset = 8 ‘Leave the y-origin 


‘Initialize x and y 
x= 8 
y=8 
’ 
‘Do 2588 iterations 
FOR n = 1 TO 2588 
pk = RND 
‘The next line works for m<=4. It must be modified 
‘for values of m> 4. 
IF pk <= p(1) THEN k = 1 ELSE IF pk <= p(2) THEN k = 2 
ELSE IF pk <= p(3) THEN k = 3 ELSE k = 4 
newx = a(k) * x + b(k) * y + e(k) 
newy = c(k) * x + d(k) * y + £(k) 
x = newx 
y = newy 
‘Use PRINT x,y instead of the PSET line 
‘to see the range of coordinates. Then fix 
‘xscale, yscale, xoffset, and yoffset 
IF n > 18 THEN PSET (x * xscale + xoffset, y * yscale 
+ yoffset) 
NEXT n 
LOCATE 24, 35 
PRINT "Press any key to end."; 
WHILE INKEY$ = "" 
WEND 
, 
‘Return to text screen 
SCREEN 8 
END 
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Figure 12: The Collage Theorem is applied to a leaf. The collage at lower left isn’t 
much good, so the corresponding IFS image, shown at lower right, is a poor 
approximation. But as the collage improves, upper left, so does the IFS image. 





Figure 13: Can you find the IFS codes for this spiral image? Only two 
transformations are needed. 





222 BYTE * JANUARY 1988 


So long as ad — cd is not 0, it is a stan- 
dard calculus result that our ratio equals 
the determinant of the transformation ma- 
trix for W,. So a good choice for the prob- 
ability p, is 


Gia — DG. 
> la,d,—b,¢,| 


provided none of these numbers p, comes 
out to be 0. A 0 value should be replaced 
by a very small positive value, such as 
0.001, and the other probabilities corre- 
spondingly adjusted to keep the sum of all 
the probabilities equal to 1. 

We now summarize the compression 
and decompression process: An input 
image is broken up into segments through 
image-processing techniques. These 
image components are looked up in the 
IFS library using the Collage Theorem, 
and their IFS codes are recorded. When 
the image is to be reconstructed, the IFS 
codes are input to the random iteration al- 
gorithm. The accuracy of the recon- 
structed image depends only on the toler- 
ance setting used during the collage 
mapping stage. 


Applications 

For graphics applications, we use a more 
sophisticated procedure that allows full- 
color images to be encoded. Combina- 
torial searching algorithms can be used to 
automate the collage mapping stage. Fig- 
ures 2, 3, and 4 were obtained using IFS 
theory at compression ratios in excess of 
10,000 to 1. These images were based on 
photographs in recent issues of National 
Geographic. A full-sequence video ani- 
mation, A Cloud Study, was shown at 
SIGGRAPH ’87. This was encoded at a 
ratio exceeding 1,000,000 to 1 and can be 
transmitted in encoded form at video 
rates over ISDN lines (ISDN stands for 
integrated services digital network, a 
concept for integrated voice and data 
communications). A frame from the ani- 
mation is shown in figure 5. 

The IFS compression technique is 
computation-intensive in both the encod- 
ing and decoding phases. Computations 
for the color images were all carried out 
on Masscomp 5600 workstations (dual 
68020-based systems) with Aurora 
graphics. Complex color images require 
about 100 hours each to encode and 30 
minutes to decode on the Masscomp. 

For practical applications, you need 
custom hardware that can speed the en- 
coding and decoding process. An experi- 
mental prototype, the IFSIS (iterated 
function system-image synthesizer), de- 
codes at the rate of several frames per 
second. The IFSIS device was produced 
from a cooperative effort between GTRC, 


IMAGE COMPRESSION 
0 EY SE RIES EPL RS et 





DARPA, Atlantic Aerospace Electronics 
Corporation, and Iterated Systems, and it 
was demonstrated on October 5, 1987, at 
the third annual meeting of the Applied and 
Computational Mathematics Program of 
DARPA. It can be connected to a personal 
computer through a serial port; the per- 
sonal computer sends the IFS codes to the 
device, which responds by producing com- 
plex color images on a monitor. 

The IFSIS is a proof of concept for 
faster devices with higher resolution. 
Once the higher-performance IFSIS de- 
vices are combined with ISDN telecom- 
munication, full-color animation at video 
rates over phone lines will be a reality. 

Another area for future application of 
IFS encoding is automatic image analy- 
sis. What’s in a picture? Does it show a 
spotted sandpiper or a robin? The more 
complex the image or the more subtle the 
question, the harder it becomes for an al- 
gorithmic answer to be formulated. But 
here’s the point: Whatever the answer, it 
will proceed faster if stable, compressed 
images are used. The reason for this is 
that image-recognition problems involve 
combinatorial searching, and searching 
times increase factorially with the size of 
the image file. 

During the spring of 1987, Iterated 
Systems was incorporated to develop 
commercial applications of IFS image 
compression. It is exciting to see how an 
abstract field of mathematics research is 
leading to new technology with implica- 
tions ranging from commercial and in- 
dustrial work to personal computing. @ 


ACKNOWLEDGMENTS 

Figures 2 through 5 were encoded by gradu- 
ate students Francois Malassenet, Laurie 
Reuter, and Arnaud Jacquin. All color 
images were produced in the Computergra- 
phical Mathematics Laboratory at Georgia 
Institute of Technology and are copyright 
1987, GTRC. 


BIBLIOGRAPHY 

Barnsley, M. F. and S. Demko. “Iterated 
Function Systems and the Global Con- 
struction of Fractals.” The Proceedings 
of the Royal Society of London, A399, 
1985, pp. 243-275. 

Barnsley, M. F., V. Ervin, D. Hardin, and 
J. Lancaster. “Solution of an Inverse 
Problem for Fractals and Other Sets.” 
Proceedings of the National Academy of 
Science, vol. 83, April 1985. 

Barnsley, M. F. Fractals Everywhere. Ac- 
ademic Press, 1988. Forthcoming. 

Elton, J. “An Ergodic Theorem for Iterated 
Maps.” Journal of Ergodic Theory and 
Dynamical Systems. Forthcoming. 

Mandelbrot, B. The Fractal Geometry of 
Nature. San Francisco, CA: W. H. Free- 
man and Co., 1982. 








Attention all FX80, FX100, JX, RX, & MX owners: 


You already own half of 





Now for $79.95 you can own the 
rest. You see, today’s new dot 
matrix printers offer a lot more 
features. 

But now, a Dots-Perfect upgrade 
kit will make your printer work like 
the new models in minutes—ata 
fraction of their cost! 

For example, with a Dots-Perfect 
the “dotty” look is gone! In its place 
is NLQ (Near Letter Quality), a 
feature that produces printed 
characters almost as sharp as a 
daisy wheel or laser printer. 
NLQ: Like Having Two Printers 


NLQis not just “double striking” 
It's a completely redesigned letter 
form that reshapes each character. 
Characters are actually formed 
by four times the normal dot 
resolution. 

And since NLQ can be switched 
on or off from the touch of a button, 
you can switch modes whenever 
you wish. 

Installs In Minutes 


You can install a Dots-Perfect in 
minutes with an ordinary screw- 
driver by following our easy-to- 
understand illustrated instructions. 

All you need do is unplug the 
existing chips and replace them 
with a Dots-Perfect—that’s all there 
is to it! 

Doesn't Get In Your Way 
ts-Perfect is easy to operate 
because it actually becomes part of 
your printer. 
It requires no software itself, so it 


can be used with every software 
package. 

And you control it from the 
printer's existing control panel 
buttons, 

You can choose from over 160 
printing modes like condensed, 
double-wide, italic, high-speed 
draft, or even combinations. All 
switched on and off at the touch 
of a button, 

Samples shown actual size 


abcdefghiJKLM 


Condensed Fine Print 


Emphas ized 


NLQ 


Doubl e-—-Wide 


Italics Underline 
Or Hundreds of Combinations! 


Dots-Perfect even has a unique 
panel button controlled buffer clear 
command. 

1.B.M. Graphics 

FX, JX and MX versions of Dots- 
Perfect will even permit switching 
from the standard Epson®character 
set to the IBM°Graphics Printer set, 
allowing you to print IBM screens 
exactly as they appear on your 
computer. All three versions have 
all Graftrax Plus™ features, even 
the MX version. 

And, every Dots-Perfect is backed 
with a full year's warranty. 

So, call now toll-free and use your 
Visa, MasterCard, or American 
Express card. 

Don't replace your printer, 
upgrade it! 


1-800-368-7737 


(Anywhere in the United States or Canada) 


=e ier 


An upgrade kit for EPSON FX, JX, RX, & MX printers is 


Circle 84 on Reader Service Card 


*&Dresselhaus 
8560 Vineyard Ave., Ste. 405, Rancho Cucamonga, CA 91730 


letter ple of 


Dots-Pertect 






1} 
(714) 945-5600 
iit IBM Is a registered trademark 
of International Business Machines Corp., 
Graftrax"™”S is a trademark of Epson America, Inc. 
son is a registered trademark of Epson America, Inc. 
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